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1. Introduction 


As an analogue of the dimension of the Riemann-Roch spaces of divisors on an 
algebraic curve, van der Geer and Schoof introduced the function for a number 
field F (see 15 ). This function is also called the “size function” for F (see 5 Hill)- 


The properties and an upper bound for /i° were provided in 15 . After that, 


10 Section 10.9], Schoof proposed a method to compute this function by using 


reduced Arakelov divisors. Essentially, we can approximate the value of at a 
given class of divisor D by knowing the short vectors of the lattice L associated 
to D because the main contributions to hP{D) come from the shortest vectors of 
L (see and [^). This can be efficiently computed if L comes with a good, i.e., 
reasonably orthogonal, basis. 

Here we present a method to approximate the value of by using the Poisson 
summation formula as well as some “good” divisors which can be obtained by 
doing the “jump algorithm” [E| Algorithm 10.8] and from an LLL-reduction on D 
(see Section]^ and 12,^). These divisors may not be reduced in the usual sense 
(see but they can be used to compute efficiently. 

Let E be a number field of degree n and discriminant Ap with unit group of 
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rank at most 2. We compute an approximate value of at any class of Arakelov 
divisors on Pici?. This method does not require a basis of the unit group of the ring 
of integers Op oi F and runs in polynomial time in log|Ap’| (see Section 5.21. 

In Section 2, we give a brief introduction to Arakelov divisors, the Arakelov 
class group and the group Picj?, the function for a number field and the Poisson 
summation formula for lattices. Then we discuss some results on Arakelov divisors 
obtained from the LLL-algorithm in Section 3. Section 4 provides an algorithm to 
approximate the values of for number fields with large discriminant. Section 5 is 
devoted to bounding the error as well determining the running time of the algorithm. 
We also present some numerical examples applying this method to compute for 
real quadratic fields and real cubic fields in Section 6. 


2. Preliminaries 

This section briefly recalls some basic definitions that are used in the next sections. 


See 10,15 for full details. 


Let F be a number field of degree n with the ring of integers Op and ri,r 2 
the number of real and complex infinite primes (or infinite places). Let Clp he the 
class group of F. Denote by Fr := F ®q M ~ Ha real complex where a runs 

through the infinite primes of F. Then Fr is an etale M-algebra with the canonical 
Euclidean structure given by the scalar product 

(M,n) := Tr(uTJ) for u = {ua-),v = (v^) € Fr. 

In particular, in terms of coordinates, we have 

||m|| 2 = Tr(uit) = ^ degcr|Mcr|^, for any u = (ua) G Fr. 


Here the degree of an inhnite prime a is equal to 1 or 2 depending on whether a is 
real or complex. 

The norm of an element u = (ucr)a of Fr is defined by N{u) := Ilcrreai'^^cr ■ 

ricT complex \'^rr\ ■ 


2.1. Arakelov divisors 

Here and in the rest of the paper we often call fractional ideals simply ‘ideals’. If 
we want to emphasize that an ideal is integral, we call it an integral ideal. 


Definition 2.1. An Arakelov divisor is a pair D = {I^u) where I is an ideal and 
u is an arbitrary unit in Yia C Fr . The Arakelov divisors of F form an additive 
group denoted by DIvr’. 


Let D = (/, u) be an Arakelov divisor. The degree of D is defined by deg{D) := 
— log {N{u)N{I)). We associate to D the lattice ul := {uf = {ucr • cr{f ))cr ■ f G 
/} C Fr with the inherited metric from Fr (see [I 10 ). For each / G /, by putting 
I1/||d := 11^/11, we obtain a scalar product on I that makes I become an ideal 
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lattice as well Section 4]. The covolume of the lattice L = ul associated to D 
is COVOl(L) = A/|Ai?| N{I)N{u) = g-degiO)^ 

To each element f G F* is attached a principal Arakelov divisor (/) = 
I/I) where f~^OF is the principal ideal generated by f~^ and |/| = 
(k(/)l)aen. R+ C Fr. It has degree 0 by the product formula. 


2.2. The group Pic^ and the Arakelov class group Pic^ 

Definition 2.2. The quotient o/Div^ by its subgroup of principal Arakelov divisors 
is denoted by PIcf- 

The set of all Arakelov divisors of degree 0 form a subgroup of Divp, denoted by 
Div^. The Arakelov class group of a number field is analogous to the Picard group 
of an algebraic curve defined as follows. 

Definition 2.3. The Arakelov class group Pic^ is the quotient ofOiv^ by its sub¬ 
group of principal divisors. 

Each V = (vcr) G 0(tK can be embedded into Div^ as the divisor Dy = {Of,u) 
with u = {e~'"‘’)a. 

Put (0£rIR)° = {(uct) G 0O-R : deg{Dy) = 0} and A = {(log |(T(e)|)o- : £ G 0*p}. 
Then A is a lattice contained in the vector space (0o.K)°. We define 

T° = (0„M)7A. 

By Dirichlet’s unit theorem, T° is a compact real torus of dimension ri 0 r 2 — 1 
Section 4.9]. The structure of Pic^ is described by the following proposition. 

Proposition 2.4. The map that sends the class of a divisor {I,u) to the class of 
the ideal I is a homomorphism from Pic^;. to the class group CIf of F. It induces 
the exact sequence 

0 — >T° ^ Pic^F —^ CIf 0. 


Proof. See 10 Proposition 2.2]. 


□ 


2.3. Metric on the Arakelov class group 


For u G Her 1^+1 we let logu denote the element logu := (log(Mcr))CT G {Icr ® 
Fr. By using the scalar product from Fr, this vector has length jjlogMjj^ = 
d-eg{a) log^(Mcr). We define 

IJullpic := min ]] logrt']] = min ]] logdEjit)]]. 

logit^=logu( mod A) 


Now let [F] and [D'] be two classes containing divisor D and D' respectively 
lying on the same connected component of Picp. Then by Proposition 2.4 there is 


some unique m G T° such that D — D' = {Of,u). We define the distance between 
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two divisor classes containing D and D' to be ||u||pic- The function || \\pic gives 
rise to a distance function that induces the natural topology of Picp. See Section 6 
in for more details. 


2.4. The function for a number field 


Let D = (/, u) be an Arakelov divisor of F. The effectivity e{D) of D is the number 
defined by 


e{D) = 



if Of ^ / 
if Op C I. 


This number is between 0 and 1. A divisor D is call effective if e{D) > 0. Similar 
to a Riemann-Roch space of an algebraic curve, we denote by H^{D) the union of 
{0} and the set of all elements f oi F for which the divisor (/) + D is effective, i.e.. 


H°{D) := {/ e F* : e((/) + O) > 0} U {0}. 


Since e((/) + O) > 0 if and only if / G / \ {0}, the set H^{D) is equal to the 
infinite group /. In order to measure its size, we weight each element / G / with 
the effectivity of the Arakelov divisor (/) + D: 

e((/) + D) = 


The value of the function hP at D is obtained by summing up these terms for all / 
in I including 0 and then taking its logarithm as follows. 


h°{D)= log 




= log 



-TrWfv 


Since two Arakelov divisors in the same class in PicF have isometric associated 
lattices, the function hP{D) only depends on the class \D] of D in PicF and we may 
write h^{[D]). In other words, is well defined on Picp. See 15 for more details. 


2.5. The Poisson summation formula for lattices 

Let L be a lattice in R" and L = Li 0 L 2 . Note that we do not require that Li and 
L 2 are orthogonal to one another. 

Consider the sum 

zeL 

Let V be the subspace Li 0 M of M" and let 7 r(b) denote the orthogonal projection 
onto F of a vector b G K". We have the following. 


Lemma 2.5. Let 7 be the covolume of the lattice Li inside V. Then 


' beL2 




E 


^ —27r2(a,7r(b)) 
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Proof. Let z G L. Then z = a + b for some a G Li and b € ^ 2 - The vectors b — 7r(b) 
and a + 7r(b) are orthogonal. The Pythagorean Theorem implies therefore 

||a + bf = ||b- 7r(b)f + ||a + 7r(b)f. 


Thus, we can write 


g — ^ ^ g—^lla+bll — ^ g—^llb—7r(b)|| ^ g-7r||a+ii-(b)|| 

h^Li2 sGLi h^Li2 aGLi 

Applying the Poisson summation formula to the second sum (see and Chapter 
7 in [^), we obtain that 


S = 


„-7r||b—n-(b)|| = 


^ —27rz(a,7r(b)) 


beL2 




□ 


3. Some Results 

In this section, we discuss “nice” properties of Arakelov divisors obtained from the 
LLL-algorithm. 

From now on, we put 

=|senK;:iV(s) = l| 

/ 2 \ "^2 _ ^ 

and dp = ( — j a/IAfI and Vp = 9p. 

Definition 3.1. Let J be an ideal of F. Then the Arakelov divisor associated to J 
is d{J) := (J, u) where u = {ua) G Ho -Ua = N{J)~^/‘^ for all a. 

Let D = (/, u) be an Arakelov divisor and let L = ul he the lattice associated 
to D. Assume that a basis of L is given. By using the LLL-algorithm, we can find 
an LLL-reduced basis {6i, ...,&„} of L. Since bi G L = ul, there is some nonzero 
element f G I such that bi = u ■ f. Denote by 

Then J is an ideal of F. Therefore, we can define as follows. 

Definition 3.2. Let D = (J, u) be an Arakelov divisor and let L = ul be the lattice 
associated to D with a known basis. We call an LLL-reduction on D the process of 
finding an LLL-reduced basis {6i,...,6„} of L, then computing a new ideal lattice 
J = bf^L and a new divisor D' = d{J). 

We first recall the following lemma Proposition 4.4]. 

Lemma 3.3. Let D = (/, u) be a divisor of degree 0. Then there is a nonzero 
element f G I such that it£r|cr(/)| < for all a. In particular \\uf\\ < ^Jn (?]/". 
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Proof. See 10 , Proposition 4.4], 


□ 


We prove the proposition below. 

Proposition 3.4. Let D = (/, u) he an Arakelov divisor of degree 0 and D' = d{J) 
obtained by an LLL-reduction on D. Then we have the following. 

i) The ideal J~^ is integral with 

Moreover, we obtain that < covol(J) < -^/j. 

ii) There is some s G and some f G I such that 

D-D' + {f) = {OF,s) 

and 


\\D — D'\\pic = ||s||pic < logPp. 

Proof, i) Since D' is obtained from an LLL-reduction on D, there is an LLL- 
reduced basis of the lattice ul associated to D such that bi = uf and 

J = /“^/ for some f G I. As f G I, the ideal J contains 1. Thus is integral. 
By Lemma |3.3[ there is a nonzero element g G J such that 

A(J)-i/"|(t( 5 )| < for all a. 

Hence 

iV(J)“^/”max|cr( 5 )| < 9^”. (3.1) 

(7 

We have 


\u{f9)\\ = \\iuf)g\\ = WhgW < max|CT(g)||16i 


Furthermore, ||6i|| < 2^” since u{fg) G ul and by the property of 

LLL-reduced bases Section 10]. As a consequence. 


max|a(5)|>^jM>2-("-i)/2. 

l|oi|| 


(3.2) 


The inequalities (3.1) and (3.2) imply that 

dp 




(max,^ \a{g)\y 




F- 


Therefore the first statement in i) is proved. Since covol(J) = ^/\Af\N{J), the 
second statement in i) follows. 

ii) The divisor D = {I,u) is of degree 0, by Lemma 3.3 there is a nonzero 
element /' in I such that ||u/'|| < ^/n 9]/". Since bi is the first vector in an LLL- 
reduced basis of the lattice ul, we again have ||&i|| < Section 10]. 

It follows that lj6il| < y/n = I?]/". 
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Let s = u\f\N{J)^^'^ = |6i|-/V(Then with the notation in i), we have that 
D — D' + (f) = {Of,s). In particular, D — D' = {Of,s) S Pic^. By Section 
this leads to the following. 


2.3 


\D-D'\ 


Pic = l|S||Pic- 


Part i) shows that J ^ is integral, so N{J) < 1. Therefore, the following in¬ 
equality holds. 


So- < ||s|| = ||6i||A^(J)^/" < 1?^" for all a. 

This leads to log(so.) < ^ logT’F for all a. Since deg(cr) logs^ = 0, we can easily 
prove the following Lemma 7.5]. 

II fogsf = X!^®s(<7)|fogs^Tp < n{n- 1) ^^logX>F^ < log^T^F- 

Since ||£> - D'\\pic = ||s||pic < || logs||, part ii) is proved. □ 


Definition 3.5. Let W = {I,v) be an Arakelov divisor of degree d. We call D = 
{I,u) with u = e'^l'^v the divisor translated from W. 

Note that if D is translated from a divisor W then deg{D) = 0. In other words, 
the class of the divisor D is in PiCp. 

We prove the following corollary. 

Corollary 3.6. Let W = (/, u) he an Arakelov divisor of degree d and let D be 
the divisor translated from W. Assume that D' = d{J) is a divisor obtained by an 
LLL-reduction on D. Then L = is the lattice associated to W 

for some s G (n<Ti^+)° iisii Pic ^ logX>F- 


Proof. By Proposition 3.4 there exists some s G (n(T®+) which ||s||Fic < 
logI?F and f G I such that D + {f) = D' + (Of,s) in Div^. Then so D = 
D' + {Of, s) = (J, iV( J)-i/”s) in Pic^. 

Since W — D + {Of, G PicF, we have the following. 

W = D + {Of, = (J, G PicF • 

Thus, the lattice associated to IF is L = □ 


4. Computing The Function h° 

Let IF = (J, v) be an Arakelov divisor of degree d and let D be the divisor translated 
from IF. Assume that a basis for the ideal lattice J and the coordinates of the vector 
V are known. We compute an approximate value of hP at the class of the Arakelov 
divisor IF G PicF with some given error 5. We consider the case in which u is a long 
vector and the discriminant Af is quite large since it is quite trivial to compute 
hP{W) in other cases. 
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We have that 

h°{W) = 

fei 

We approximate the value of h^{W) with some small error. This can be done 
by summing up only the large terms, i.e., the terms for which ||/f |P < M 

with some given M > 0. In case n is a long vector, while collecting such short 
vectors /, it is quite easy to miss many of them. Consequently, the obtained value 
of h^{W) may be smaller than the true value. Therefore, we will find some “good” 
divisor D\ that is obtained from an LLL-reduction on D and has nice properties 
described in Section]^ then use it for computing h^{W). 

Note that by Proposition |3.4[ for any given divisor D of degree 0, there exists a 
good divisor D' close to D in Pic^;. in the sense that \\D — D'Wpi^ < logVp. 

We first describle the following algorithm that is similar to (cf.[ [^, Algorithm 
10.4]). 


ALGORITHM 4.1. Given two Arakelov divisors Di = d{Ji) and D 2 = d{J 2 ) such 
that iV(Jfand N{J^^) < , compute a divisor d(J) 

obtained from the LLL reduction on D\ +D 2 S Picp in polynomial time in log |Ai?|. 
Description. Since N[J^^) < and the result 

D 3 = Di + D 2 = {J 1 J 2 , N{JiJ 2 )) can be computed in time polynomial in log |Ai?|. 
Then one performs the LLL reduction on the divisor ZI 3 . The resulting divisor d{J) 
is then close to Di + D 2 by Proposition 3.4 Since N{JiJ 2 )~^ < 2'^^^~^'>dp, the 


_ Since ^(JiJa)-! 

running time of this second step is also polynomial in log jAi^’j 


Next, we explain how to compute efficiently a divisor D' obtained from some LLL 
reduction close to a given divisor D = {Of,u) in Pic^. This process can be seen as 
performing repeatedly doubling and LLL-reduction to go from the origin (Of, 1) to 
D. We apply the “jump algorithm” [10, Algorithm 10.8] with a minor modification 
to adapt to our situation. Indeed, in the reduction step, instead of using a shortest 
vector, we use the first vector of an LLL-reduced basis of the lattices associated to 
Arakelov divisors. 


ALGORITHM 4.2. Given a divisor D = {Of,u) of degree 0, compute a reduced 
Arakelov divisor whose image in Pic^ has distance less than logI?F from D. 
Description. Assume that u = (e““'‘')cr- Let t > 0 be the smallest integer for which 
n • 2“‘ • \wa\ < logc^F for all cr. Then = 2“* • Wa satisfies n ■ \za\ < log^F for all 
cr. Let u! = (e~^‘')a-. Then = u'. In other words, from the point (Of,w) we can 
reach to D after t times doubling. Denote by Wi = (Of,uj‘^ ) for i = 1,2, ...,t. We 
inductively compute Arakelov divisors D[ = d{Ji) obtained by LLL reduction for 
which 

\m-D[\\p,,<\ogVF. (4.1) 

We compute from Wi by doubling and doing LLL reduction. More precisely, 
by induction, there exists some uJt G such that Wi = D'^ + {OF,^i) in 
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Pic^;^ and 


\Wi ~ O'iWpic = llWillpic < logPp. 


Since Wi = D'^ + {Of,uj,) = {J^,Vi) where Vi = we get W,+i = 

2Wi = {Jf, vf). Let be a divisor obtained by an LLL reduction on Wi+i. Then 
there is an LLL-reduced basis of the lattice Li+i = vfjf associated 

to Wi+i such that J^+i = Proposition 

(OpjWi+i) in Pic^ for some G (nCT®+)° 


shows that Wi+i = + 


||fPz+l — D'^+i\\pic = ||Wi+i||pic < logPp. 

Note that here uji+i = |6i|-/V( Thus, we can construct all divisor D[ satis¬ 

fying (4.11 for j = 1, 2,..., t. 

Now let s = (jJt and let D' = d[Jt). Then 


D = Wt = D'^ + {OF,u:t)=D' + (Op, s) in Pic^. 

for s = ujt & (11,7®+)° and \\D - D'\\pi^ = ||s||pjc = \\i^t\\pic < logT>p. This 
completes the description of the algorithm. 

With the notation of Proposition |3.4| and Corollary |3.6[ let L = 
e-d/"-]V(j)-i/’T’gj be the lattice associated to W. Then 

h°(W) = log^ = log 5] (4.2) 

fGJ 9'^^ 


Every vector of the lattice L has the form e ‘^/”A^(J) for some f € J. 

As s is short, 1V(J)“^/" is a small scalar and J is a nice lattice (see Proposition 


the lattice L more efficiently. Therefore, the value of h'^{W) can be computed more 


in succession. 


exactly. Computing h^{W) is done in 3 steps described in Section 


4.1 


4.2 


and 4.3 


3.4), we can easily compute an LLL-reduced basis of L then find short vectors of 


4.1. Finding a good divisor D' close to D 


Assume that a basis for the ideal lattice I and the coordinates of the vector v are 
known. We will find a divisor D' that is obtained from some LLL reduction and 
with the property that D — D' + (Op,s) in Picp for some s G (11^®+) and 
p-^'||«c<31ogI?p. 

Let Di = and to £>2 = (Op, N{iy/^u). Then Di, D 2 have degree 

0 and D = D 1 +D 2 . We compute divisors D[ (see|4.1.l]) and (see 4.1.2) obtained 


from some LLL reduction so that ||Zli — Zl'i||pic < logUp and ||Il 2 — < 


log (Dp. Then we find a divisor D' = d{J) close D'^ + D 2 in PiCp (see 4.1.3). This 
process is described as follows. 
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4.1.1. Computing Hj. 

We compute a divisor D'^ close to Di in Pic^^ in the sense that its distance to Di 
is at most logP p. T his can be done easily by performing the LLL reduction on Di. 


By Proposition 


Picp and 


3.4 


there is some si S {U.K) so that Di — D[ = {Of,si) in 


\Di — D[\\pii. — ||si||pic < ^ogVp. 


4.1.2. Computing D 2 

Use Algorithm 4.2 to compute a reduced Arakelov divisorU^ whose image in Pic^;. 


has distance less than logVp from D 2 = [Of,u') with u' = N{iy/'^u. We obtain 
that 

£>2 = £>2 + (Op, S 2 ) in Pic^^ 


for S2 e (n,7®+)° and \\D2 - D'^Wpic = ||s2||pjc < logX>p. 


4.1.3. Computing D' 

Adding divisors D'^ and £>2 as described in Algorithm |4.1[ we then compute a divisor 
D' — d{J) close to D'l + £>2 in PiCp. Indeed, by performing LLL reduction on the 
divisor D'^ + £>2 we obtain D' = d{J) and {D[ + D' 2 ) — D' = (Op, S 3 ) for some 
■S 3 e (11^®+)° and \\{D[ + D' 2 ) - H'ljpic = ||s 3 ||pic < logX>p. 

Let s = Si • S 2 • S 3 G (no-®+)°- Then 

D-D' = (D 3 - D'^) + (D 2 - D' 2 ) - ((D'^ + D' 2 ) - D') = (Op, s) G Pic^ . 

Thus, \\D - D'Wpic = ||s||pjc < ||si||pic + \\s 2 \\pic + llssllpic < 31og£>p. 


4.2. Applying Poisson summation for the lattice L 

Since D — D' = (Op,s) G Pic^, it follows that D = (J, A^( J)“^/"s). By translating 
D to the divisor W, we obtain that W = D + (Op, = (J, J)“^/"s) . 

Let £ = J)“^/”sJ be the lattice associated to D. Assume that £ has an 

LLL-reduced basis {bi,..., b„}. Let {b*,..., b* } denote the Gram-Schmidt orthogo- 
nalization of this basis and 


Pi,j — 


(b^b*) 

l|b;p 


for all 2 < i < n and 1 < j < i — I. 


Then any element z G L can be written uniquely as z = with the 

coefficients Xi G Z for all i = 1, 2,..., n. Similar to Lemma [4^ we can write |jz|p as 
below. 

n / i— 1 

l|z|l^ = I + 

i=l \ j = l 



q(Xi,X2, .■.,Xn). 


(4.3) 
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where = |lb*|p and Aij = Therefore, 

/lO(IF) = (4.4) 

z£L XiGlj 

Remark 4.3. We only catch the vectors z in L for which ||z|p < M. In other 
words, we only compute vectors x = {xi,X2, ■■■,Xn) S Z" satisfying q{x) < M. The 
Fincke-Pohst algorithm Algorithm 2.12] or an LLL reduced basis of L Section 
12 ] can be used to find the list C of these vectors x. 

An approximate value of is obtained by summing up only the terms 
for which x S £ as below. 

/i°(IF) « log^e-’^'jW. (4.5) 

xGC 

The lattice L has covolume covol(L) = Jir^i = VTli=i ^ 

have at most a(n) vectors. Here a{n) is a function depending only on n. See 

Algorithm 2.12 in and Section 12 in for the explanation. Therefore, in order 

to reduce the number of vectors in the list £, we can “make” covol(£) larger by 
using the Poisson summation formula as follows. 

Assume that = jjbijj < 1. Let k be the largest index such that Ai^i < 
1 for alH < k. Denote by 

Ai = ©f=iZ • bi and £2 = ©"=fc+iZ • bj. 

Then L = £1 © £ 2 . Since {bi,..., b„} is LLL-reduced, the vectors bi,..., bk form an 
LLL-reduced basis for £1 (see [^). 


Remark 4.4. Let B be the matrix of which columns are vectors bi,..., b^ and let 
G = B* ■ B. Then the columns of B ■ G~^ form a basis for the dual lattice L'( of £ 1 . 


Now we apply the Poisson summation formula for £. See Lemma 2.5 Let 7 be 


the covolume of the lattice £1 inside F = £ 1 ' 


L Then 7 = n*=i l|b -j] and 


h°{W) = log I ;J; XI X 

beL2 aGLj' 


^ —27rz(a,7r(b)) 


(4.6) 


Assume that £]^ has a basis Ci, • • • , that is computed by Remark 4.4 Denote 
by {c]], • • • , c^} the Gram-Schmidt orthogonalization of the basis {ci, • • • , Cfc} and 


a,; = lie*! 


= 


(c*,cp 

llr*l|2 


for all 1 < i < fc and 1 < j < t — 1 . 


Now let a = ^ £1 where Xi G Z for all i = 1,2, ••• ,fc and b = 

Sj=fc+i 2^1 bj G £2 where xj G Z for all j = A: + 1, • • • , n. 

Lemma 4.5. We have 




for all 1 < i < k 
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and moreover 


= Ci^i ( + X! ‘^’■4 


(4.7) 




Proof. This is easily proved by using Remark 4.4 and properties of the Gram- 
Schmidt orthogonal basis {c);, • • • , c^}. □ 

Recall that Ajj = ||b*|p and Atj = fj-tj for all j = fc + 1, • • • ,n and t > j. We 
have the lemmas below. 

Lemma 4.6. We have 


k n k 

(a,7r(b))=^ ^Aj-i(Q,b*)a:/Xj. 

j^k-\-l i^l 


(4.8) 


Proof. Because 7r(b) is the orthogonal projection of b on F that has an orthogonal 
basis b);, • • • , b^., we obtain that 


k ,, , k n 

= E = E E 

i=l i=l j=k+l 


(4.9) 


Then the result is implied by taking scalar product of 7r(b) with a = 
Lemma 4.7. We have 


|b-7r(b)f = Y ^ 




j=k+l 


t=j+l 


Proof. Since bj = b* + a.11 j > 2, the vector b therefore can be 

rewritten as 


E + E fitjXt b* + E E 


(4.10) 


j=fc+i 


t=j+i 


i=l j—k-\-l 


By using equalities (4.71 and (4.10), the result is obtained since the vectors 
b^_,_^,..., b* are pairwise orthogonal. □ 


Lemma 4.5 4.6 4.7 and (4.6) lead to 


hO(W)=log -^e-^ 


—'kQ{Xi,X 2 Xji) 




(4.11) 
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where Q{xi, ■ ■ ■ ,a;„) = Qi(xi,x„) + 2Q2(a;i, • • • ,Xn)i with 
Qi{xi,...,Xn) = ||b - 7r(b)f + |ja|p 




k / k 

j + ^ Ci.i (+ X! Cr,iXr ) (4.12) 

t=i+i 


i=l 


r—i-\-l 


and 


k n k 

Q 2 (a:i,-• • ,x„) = (a,7r(b)) = ^ ^ '^Aj^i{ci,h*)xiXj. 

1=1 j=k+l i=l 

4.3. Finding the short vectors of the lattice associated to D 

An approximation of h'^{W) is obtained by summing up the terms such 

that Qi(x) < M. By using the Fincke-Pohst algorithm, we can find the list £i 
containing all vectors x = (cci, • • • ,a;„) € Z” such that Qi(x) < M. See Algorithm 
2.12 in 1^. Then an approximate value of hP{W) is obtained as follows. 


/i“(IF) « log 




(4.13) 


xGz;i 


Remark 4.8. Let £i = {x € Z" : Qi(x) < M} and £ = {x G Z" : g(x) < M}. Let 
L' be the lattice associated to the quadratic form Qi(x). Then 


covol(£')^ = n ^*4 II 

2=1 j—k-\-l 


- -pffe . n 

ll2 = l j = fc+l 


, ■ covol(£)^. 

(ntiAO^ 

(4.14) 

We have that #£ < ciin-) and #£i < CK{n) (see Algorithm 2.12 

in and Section 12 in !)• Since Ai^i < 1 for all j < fc, it follows that covol(£') > 
covol(£). From this inequality and (4.141, we usually obtain that #£i < #£. 

Thus, the list £i usually contains less vectors than the list £ does. In addition, 
1 _ 1 


Vn, 


> 1. At the result, the sum in (4.13) converges better than in (4.5). 


Hence, we can compute (4.13) by only summing a small number of terms. 


Note that the function qfminim in pari-gp that uses the Fincke-Pohst algo¬ 
rithm, can be used to find all nonzero vectors (up to a sign) with length bounded 
by M of a given lattice. Another method uses an LLL reduced basis of the lattice 
L'] see Section 12 in [^. For a fixed lattice, the complexity of both methods is in 
polynomial time in M (see Section 5.2). 

By the proofs of 

to approximate hP{W) with an error (5, we 


5.1 


5.2 


5.3 


and Proposition 

(log(l/'^) + (’^ + 1) log 3 -h (n(n + l)/2- 1) log 2). 


Lemma 
can choose M 

The algorithm below computes an approximate value of h^{W) with a given 
error <5. 

Input: 
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• A basis for the lattice I. 

• The coordinates of v. 

• An error 5. 


Output: An approximate value of h^{W) with error 5. 


4.1 


ALGORITHM 4.9. 

(1) Find a divisor D' that is close to D in Pic^^- as described in Section 

(2) Apply Poisson summation formula. 

(a) Find an LLL-reduced basis {bi,..., b„} of L'. 

(b) Compute {b*,...,b*} and Ai^i = ||b*|p and Aij = for ah 2 < 

i < n and 1 < j < i — 1. 

(c) If ||bi|| > I, then put Qixi, ...,a:„) = Qi{xi, ...,Xn) = qixi, ...,x„) (see 


(4.3)) and <32(2:1, L' = L and 7 = 1. 

If IIbill < 1, then let k be the largest index such that ||b*|| < 
1 for all j < k. Denote by B the matrix of which columns are vectors 

bi,..., bfc. 

(i) Compute G = B* ■ B and C = B ■ G~^. 

(ii) Let Cl, • • • , c/j be the columns of G. Compute {cj, • • • , c^} and Gi^i = 

||c*|| and Gij = for all 1 < i < fc and 1 < j < i — 1- 

(iii) Compute (c;, b*) for all I = I,..., k and f = 1,..., k. 


(iv) Denote Q{xi, ..., x„), Qi{xi, ..., x„) and Q 2 {xi, ..., x„) as in (4.12) and 
let L' be the lattice associated to <3i(x) and 7 = yTliLi^^^. 

(3) Find the short vectors of the lattice L'. 

(a) Compute M = (log(l/(5) + (n + 1) log 3 + (n(n + l)/2 — 1) log 2). 

(b) Find the list £1 = {x = (xi, • • • ,x„) S Z" : (3i(x) < M} and approxi¬ 


mate h'^{W) as (4.13). 


5. The Error and Running Time of The Algorithm 
5.1. Bound for the error in Algorithm \4-S\ 

To find a bound for the error in approximating the value of hP{W) in Algorithm 


Lemma 5.1. Let L' be a lattice of rank n. Assume that the length of shortest vector 
of the lattice L' is A and M > max{A^, ^ log ^}. 

Let 


4.9 we use the idea of Section 4] as below. 


S'= E 


e 


aGL' 

||a|p>M 


2 


-TT a| 
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Then In particular, the hound for goes to zero when M 

tends to infinity. 

Proof. Let Bt = {s G L' : M < ||a|p < t} for each t > M. The balls with centers 
di G Bt and radius A/2 are disjoint. Their union is contained in the (hyper) annular 
disk 

{z G Fr : V M — A/2 < ||z|| < ^/t + A/2}. 


Consequently, the following is implied. 



The second inequality is since t > M > A^. Using this inequality, we get 


s'= U 


3^1/ 

||a|P>M 


7re dt < TT / fiBtC dt 

^ J M 


< TT 


IM 




! M 


( 2VM 


- 1 


dt. 


Since M > | log we have < (|)" e*. This implies that the first integral 

is at most (|)" ^ xhe second one is equal to ^ — l) e“ 


^ — 7tM 


Hence 


5" < 




/ 2y/M 


- 1 


— -nM 


TT — 1 \ A , 

Thus, the lemma is proved. □ 

Lemma 5.2. Let L' be the lattice otK" associated to the definite positive quadratic 
form (5i(x) in (4.12) of Section Then the shortest vector of the lattice L' has 
length A > 2(“"+i)/2^ 


Proof. 

If Hi 4 = ||6i|p > 1, then by the property of LLL-reduced bases, the length of 
the shortest vector of L' is at least |j Section 10]. In other words, 

A > 

If Hi i < 1, then > 1 for alH < fc and > I since k is the largest 

index such that < 1 for all i < k. On the other hand, we have = |jb*|p for 
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all i = l,2,...,n. So, if fc + 2 < j < n then Ajj = ||b*|p > = 

Ak+i^k+i > 2“"’+^ Section 10]. Thus, all the coefficients Ci^i = 
with i <k and Ajj with j > fc + 1 are at least 2“"+^. As the result, 

Qi(x) > min{Ci^i, Ajj : 1 < i < k and k + l<j<n}> 2 “"”''^. 

The result now follows since = min{Qi(x) : x G Z}. □ 


Proposition 5.3. Let S be the error in approximating hP{W) described in Algo¬ 
rithm 4-9 Then for fix degree of the number field, and for M > niax{A^, ^ log ^}, 


we have 6 = 0(e 


Proof. Now let 


E 


„—rrQ(x) 


and 


^0= E 


,-irQ(x) 


xGZ" 

Qi(x)<M 


Then h^{W) = log(i5'o) and we approximate it by log(iS'). The error in approxi¬ 
mating hP{W) is 


5 = 


log ( -5o 


log ( -S 


log 


Furthermore, 


l^o-^l = 


E 


„-7rQ(x) 


xGZ"- 

Ql(x)>M 


^ E 


„-7I-Q(x)| _ 


xGZ" 

Qi(x)>M 


= E - 

xGZ”' 

Qi(x)>M 


„-7rQi(x) _ 


= E 


sG L 

||a||2>M 


Since M > max{A^, ^ log |} > 2^^ "1, Lemma 
S' = 

Since 5'o > -S' > 1, it follows that 

So 


5.1 


and 


5.2 


show that jS'o — 5] < 


(5 = 


log 


S' 


< |S-So| < 0(6-1’^-^!^'^). 


□ 


5.2. Run time of Algorithm 

We prove the proposition below. 

Prop 5.4. Let 0 < (5 < 1. Assume that the given basis of the ideal lattice / and 
the vector v have size bounded by If the degree n of the number field is 


fixed, then Algorithm 4.9 with an error S runs in time in poly(log(l/(5) • log lAp-j). 


Proof. The basis {bi, ...,5„} of the ideal lattice I and vector v have size at most 
Therefore, Step l-finding a good divisor D' close to D by using the “jump 


algorithm”-runs in time polynomial in logjAi?] 10 Algorithm 10.8]. In addition, 
the entries of the matrix B bounded by jAi?j‘^lb since they are coordinates of 
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{bi,...,bn}- Thus, each Step 2a), 2b) and 2c) and hence Step 2 can be done in 
polynomial time in log |Ai?|. 

The list = {x = {xi,-- ■ ,x„) S Z" : Qi(x) < M} can be computed by the 
Fincke-Pohst algorithm (see Algorithm 2.12 in [^). If the degree n of the number 
field is fixed, then the complexity of this algorithm is at most 1/ covol(L') •poly(M). 
See Section 3 in for more details. The covolume of the lattice L' is covol(L') = 
ni=i,fc ■ T\j=k+i, ^ ||b*||. By a similar argument in the proof of Lemma |5.2l 
one can show that ||c*|| = Ci^i > 1 and ||b*|| = for all 

I < i < k and k + 1 < j < n. Consequently, covol(L') > Therefore, 

the complexity of Step 3 is bounded by poly(M). Proposition |5.3| says that M is 
bounded by 0(log(l/J)). As the result. Step 3 can be done in time poly(log(l/(5)). 

Overall, the algorithm runs in time in poly(log( 1/(5) • log jAi?!) for fixed degree 
n of the number held. □ 


6. Some Numerical Examples 

We compute the value of for real quadratic helds and number helds with unit 
group of rank 2. In the examples below, we pick an irreducible polynomial P of large 
discriminant and compute for the number held F dehned by P. The algorithm 
works well without requiring the units of F. Here pari - gp is used to compute 
approximate values of and Mathematica is used to plot it. 

Since the symmetry induced by Riemann-Roch (see Proposition 1 in 15 ), the 
graphs of on the cosets of PicF are similar. See Example 1, 2 and 3 in 15 . In the 
following examples, we compute on the coset Pic^^ of class of divisors of degree 
d = l/ 21 og Aj’. 


Example 6.1. Let Ap = 10®° + 129 and P = — A^? be the polynomial dehning 

F = Q{^/Ap). Then F is a real quadratic held with the discriminant Ap and with 
two real inhnite primes CTi : F —> sends \/ Ap to itself and (T2 : F —y R^ sends 
Ap to — Ap. 

The class number of F is 1 and the group Picj’ is isomorphic to a cylinder. Eor 
every d G R, the coset Pic^^ of classes of degree d is a circle whose circumference 
is equal to the regulator Rp of F. 

We have (R^) = {{xi,X 2 ) : xi,X 2 G R and xi +X 2 = 0} = {(—xi, xi) : xi G R} 
is the bisector of the second quadrant of the axes. It is a 1-dimensional subspace of 
R^ with an orthonormal basis e = (—;^)- The connected component of identity 

of Pic^ is a circle T° = (R^)°/A where A is the lattice {(log |cri(/)|, log |cr 2 (/)|) : 
/ G Op} (that is unknown). 

Denote hy w = 10^° • e = = iwi,W 2 ). Let W = {Op,v) G Pic^^ 

where v = \Ap\~^/^ ■ exp(r(;) = lA^I^^A . (e““be“’"^). The divisor D translated 
from W is {Op,u) G T° with u = • v = (e“™be““'^) and ||r(;|| = 10^° can be 

seen as the distance from iOp^ 1) to D. 

Input: 











March 23, 2016 4:21 WSPC/INSTRUCTION FILE thetafunction 


18 Ha Thanh Nguyen Tran 

• A basis for the lattice Op-. ai = (1,1) and 32 = (—5^® — 0.5,5^® — 0.5). 

. V = |A;.|-i/4 . (e-“Se-“=) and = (-^, ^). 

• 10 -^ 


We apply Algorithm |4.9| as follows. 


( 2 ) 


(1) Find a divisor D' close to D in Pic^^. 

The finite part of IF \s I = Op. As the notations in Section [TT] Di is the zero 
divisor {Op, 1), D 2 = D and u' = u. We can skip part 4.1.1 and part 4.1.3 in 
Section l4.1l To find a divisor D' obtained from some LLL reduction close to 

do Algorithm |4.2[ as follows. 


the it is sufficient to do part 4.1.2 


i.e., 


The smallest integer t such that n • 2 * • < logdi:’ for i = 1, 2 is f = 61. 

Let {zi,Z 2 ) = 2-61 .re « (-30.66587,30.66587) and w = (e-^Le-^^). 


As described in Algorithm 4.2 denote by Wi = {Op,uj'^ ) and Wi = D{ + 


{Op,LOi) with D{ = d{Ji) a good divisor obtained from an LLL-reduction on 
Wi- By performing doubling and LLL-reduction 61 times, we can reach to D. 
The result can be seen in Table [H 

In Table the second column contains the matrices W for which Mi = 
Mq ■ Ni where Mq and Mi are the matrices of which columns form a basis of 
Op and Ji respectively for all i = 1,2, ..., 61. 

Let J = Jei and s = wgi « (g-o.80975^ g0.80975^_ Then by choosing D' = 
= d{J), we obtain that D = D' + {Op, s) = {J, sN{J)~^^") G Pic^ . 
Apply Poisson summation formula. 

Since D = (J, sA( J)-i/") e Pic^ and • A(J)-i/" = (covol( J))-!/^ « 

13.46966, we obtain that IF = I? -I- {Op,e~‘^^'^) = (J, s(covol( J))”^/^). The 
lattice associated to IF is L = (covol( J))-i/"sJ. 


0.38966 and A2 2 = 


(a) L has an LLL-reduced basis {bi,b 2 } with 
61 « (0.12124,0.61234), 62 « (-1.57394,0.29870). 

(b) 6* = 61, b*2 = (-1.57148,0.31114), Ai,i = ||b)[||^ 

||b^||2 « 2.56635, A2.1 = {h2,b*i)/\\blf « -0.02033. 

(c) Since ||bi|| < 1 and Hb^H > 1, w e have k = 1. Let Li = Z-bi and L 2 = Z-b 2 

2.5 Then B = (5i) = 


as the notations in Section 


0.61234/■ 


(i) G = B^ ■B={ 0.38966 and C = B • G-^ = 


'0.31115\ 

1.57147/■ 

(ii) Cl = (0.31115,1.57147) is the column of C and c* = (0.31115,1.57147), 
Gi,i = llctip = 2.56635. 

(iii) (ci, b*) = 1 and (ci, bj) = 0. 

(iv) Qi{xi,X 2 ) = 2.56635a;i + 2.56635a;2, <32(2:1,X2) = —0.02033a:ia:2 and 
Q{xi,X2) = Qi{xi,X2) + 2Q2{xi,X2)i. 

Let L' be the lattice associated to (3i(x). 

(3) Find short vectors of the lattice L'. 
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(a) To approximate h^{W) for quadratic fields with an error 5 ~ 10“®, it is 
sufficient to choose M = 8. 

(b) The Fincke-Pohst method (Algorithm 2.12 in [^) is used to find the list 
£2 of all columns vectors x = {xi,X 2 ) G Z^\{0,0} such that Qi(x) = 
2.56635a;f + 2.56635x2 — ^■ N^ote that the function qfminim in pari-gp 
can be used to hnd £2- Here £2 contains only 4 vectors (up to a sign). 


£2 


0 111 

10 1-1 


By symmetry and since (5(0,1) = (5(1,0), an approximate value of h^{W) 
is obtained as follows. 


h°{W) 


VVO.38966 V JJ 


= 0.47250. 

Output: h^{W) « 0.47250 with an error S = 10“^. 

Recall that the coset Pic^^ is a circle containing the point W. Let X = W and 
let y = W + {Of,v') be the points on Pic^^ where log(w') = 50 • e. To see what 
hP looks like, we compute at the points in the interval [A, Y\ on Pic^ and plot 
it. The points in this interval have corresponding translated divisors in the interval 
[A, B] on where A = {Op, u) and B = A + {Op, v'). 

First, the interval [A, R] is divided into small intervals of length 1. After that, 
we do LLL-reduction on the middle points of the small intervals to obtain good 
divisors. Let S be the set of all good divisors obtained by this way. Then S has 18 
divisors in total (see Figure [^. 



Now, let W be an arbitrary divisor in [A, Y], Then its translated divisor I? is a 
divisor in [A, B]. We search for a good divisor D' in S which is the closest to D and 
use D' to compute hP{W'). Then is plotted in the interval [A, Y] as in Figures]^ 
below. In which the red points are divisors on Pic^^ whose translated divisors are 
divisors in S. 
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A ; . A 


* • •• A • . . 

• ■ : * * A A / '• • A. / \; : / ■ ; *. a 

i W *V* V Vv V V* V V V* \Af* 


Fig. 2: for the quadratic field (5(^10®*^ + 129). 


Example 6.2. Let P = - 88998X2 _ i090173446X - 1000470997815. Then 

P is an irreducible polynomial with 3 real roots denoted by ai = a,a2;Q^3- 
Let F = Q(a). Thus, F is a real cubic field with 3 real infinite primes Ui : 
F —sending a to for i = 1,2,3. The discriminant of F is Aj;’ = 
10000820940380105429207549453 > lO^s. 

We have (M^) = {(a:i, a;2, X3) S : a:i+a;2 + a;3 = 0} a 2-dimensional subspace 
of with an orthonormal basis {ei = (^.0,-^),e2 = 

The connected component of identity of Pic^;^ is the Dirichlet torus T° = (K^)° / 
A where A is the lattice {(log |cti(/)|, log |cr2(/)|, log |cr3(/)|) : / e 0*p}. 

Denote by w = 10^° • ei + 10^° • 62 

« (11153550716.50411,-8164965809.27726, -2988584907.22685) = (wi,W2,W3). 
Let W € where v = • exp(w) = IAfT^/® • (e-“b e-“’^ The 

divisor D translated from W is {Of,u) G T° with u = -v = {e~^^ ,e~^^ 
and ||ic|| « 14142135623.73095 can be considered as the distance from {Op,!) to 
D. 


Input: 


• A basis for the lattice Op-. ai = (1,1,1), 

32 « (-39667.47830,-30666.32509,70333.80339), 

33 « (451263681.71217,-491488332.71468,40224652.00251). 

• V = |Af|“^^® • (e““b6““=*) where 

{wi,W 2 , W 3 ) « (11153550716.50411, -8164965809.27726, -2988584907.22685). 

• 10 -^ 


We compute h^{W) by using Algorithm 


4.9 
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( 1 ) 


Find a divisor D' close to D in Picp 


As the notations in Section 4.1 Di = (OpA) the zero divisor, D 2 = D and 


u' = u. Similar to Example 6.1 we can skip part 4.1.1 and part 4.1.3 To find 
a divisor D' obtained from some LLL reduction close to the D, it is sufficient 
to do part |4.1.2| i.e., do Algorithm |4. 2 [ as follows. 

We have that t = 30 is the smallest integer for which n ■ 2“* • |tCi| < logdi?’ 
for i = 1,2,3. Let (21,02,-23) = 2-^° ■ w ^ (10.38755,-7.60422,-2.78333) and 


= {^ 


-21 ^-22 ^-23 


)• 


As described in Algorithm 4.2 denote by Wi = (Of,uj^ ) and Wi = O- + 


(Op,uJi) with D[ = d{Ji), a good divisor obtained from an LLL-reduction on 

A. 


In Table the second column contains the matrices of which columns form 
an LLL-reduced basis for the lattices for all* = 1, 2, ..., 30. The coordinates 
of these vectors are computed with respect to the basis of Of- 

Let J = J 30 and s = wso « (gi-09384^ g-o.6i605^ g-o.47780^_ Then by choosing 
D' = Ogg = d( J), we obtain that D — D' + {Op, s) = (J, sA( J)“^/") G PIcf ■ 
(2) Apply Poisson summation formula. 

Since D = (J, sA(J)-i/") G Pic^ and = (covol(J))-!/^, it 

follows that W = D + {Op,e~‘^^^) = (J, s(covol(J))“^/^). Here = 

938139713086 and covol(J) « 106.59831. The lattice associated to IT is L = 
(covol( J))“^/"sJ. 


(a) L has an LLL-reduced basis {bi, 62, 63} with 

bi « (0.07064,0.39051,0.34009), b2 « (0.83795,0.53263,-0.87506), 
ba « (1.35900,-0.65069,0.22025). 

(b) 5* = 61, 5* = (0.84581,0.57611,-0.83720) and 
bl = (1.09497,-0.72624,0.60648). 

Ai,i = ||b^||2 « 0.27314, A2.2 = ||b5||2 « 1.74820, A3,3 = Ijb^f « 2.09420 
and A2.1 = {h 2 ,bl)/\\btf « -0.11133, Aa.i = {b 3 ,bl)/\\blf « -0.30460 
and A3,2 = (b3,6^)/||62f« 0.33761. 

Let Li = Z • bi and L2 = ^ • b2 0 Z • 63. 

/ 0.07064 \ 

(c) Since |jbi|| < 1 and Hb^H > 1, we have k = l and B — (^bi) = j 0.39051 j . 

V 0.34009/ 


/ 0.258609 \ 

(i) G = B* 3 = (0.27314) and C = B ■ G-^ = 1.42968 . 

V 1.24508 / 

(ii) Cl = (0.258609,1.42968,1.24508) is the column of C and 
c* = (0.258609,1.42968,1.24508), Ci,! = ||c*f = 3.66108. 

(in) (ci, bj) = 1 and (ci, b^) = (ci, bj) = 0. 

(iv) Denote by Qi{xi,X 2 ,X 3 ) = 3.66108x1 + 1.74820(x2 + 

0 . 33761 x 3 )^+ 2.09420x|, Q 2 {xi,X 2 ,X 3 ) = 0.11133xiX2+ 0.30460xiX3 and 

Q(xi,X2,X3) = Qi(xi,X2,X3) 0 2(22 (a^l , 2^2 , X3)L 
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Let L’ be the lattice associated to Qi(x). 

(3) Find short vectors of the lattice L'. 


(a) By Proposition 5.3 to approximate h^{W) for cubic fields with an error 
5 ~ 10“^, it is sufficient to choose M = 9. 

(b) The Fincke-Pohst method (Algorithm 2.12 in 0 ) is used to find the list 
Li of all column vectors x = {x\,X 2 ,Xz) € Z^\{0,0, 0} such that Qi(x) < 
M. Here Li has 17 vectors (up to a sign) obtained by using the function 
qfminim in pari-gp. 

/OOOOO 0 00111111 1 1 1 

£2=1201-1-2-1001-101-1 0 1 -1 
VOOll 1 1 2 2 0 0 0 1 1 1 -1 -1 - 1 , 

By symmetry, an approximate value of hP{W) is obtained as follows. 


h°{W) « log 


1 


V0.27314 


I + 2E 


7rQ(x) 


= 0.65882. 


xG>C2 


Output: hP{W) « 0.65882 with error (5 = 10“^. 

Putting vi = 1, vi = exp (6 • ei), V 3 = exp (6 • ei + 8 • 62 ) and V 4 = exp (8 • 62 ). 
Let Xi = W + (Of, Vi) be the points on Pic^^ for i = 1, 2, 3,4. Similar to Example 
6.1 to see what looks like, we compute at the points in the box X1X1X3X4 
on Pic^ and plot it. The points in this box have corresponding translated divisors in 
the rectangle ABCE on r°. Here A = {Op, u),B = A + {Op, vi),C = A + {Op, v^) 
and E = A + {Op, V4) . 


E C 



Fig. 3: The set S of good divisors in the rectangle ABCE. 


The rectangle ABCE is divided into small squares, each one has sides of length 
1. After that, we perform LLL-reduction at the center of such squares to obtain 
good divisors. Let S be the set of all good divisors obtained by this way. Then S 
has 15 points in total (see Figure [^. 
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Now, let W be an arbitrary divisor in X 1 X 2 X 3 X 4 . Then its translated divisor 
H is a divisor in ABCE. We search for a good divisor D' in S which is the closest 
to D and use D' to compute h?{W'). Then /i° is plotted as in Figure]^ andbelow 
in which the red points are divisors on Pic^^ whose translated divisors are good 
divisors in S. In Figure]^ the dark color area is corespondent to the large values of 
h°. 




Fig. 5: /i° for the real cubic field F. 
















Table 1 


to 

4 ^ 



742409068975056334669660076059992618001 (0.80975, -0.80975) 


March 23, 2016 4:21 WSPC/INSTRUCTION FILE thetafunction 



Table 2 


i 

Ji 

N{Ji logWi 

0 

Of 

1 (10.38755,-7.60422,-2.78334) 

1 1 1 

/ 691582920399 118502077239 -227460374946^ 

0 40389624 42667155 

V 0 3277 -13661 y 

1 691582920399 (1.05601,-1.5469,0.49083) 

2 ^ 

/ 222208932162 82987136358 68645452464 \ 

r 0 19963734 42669156 

\ 0 4774 -927 / 

222208932162 (-0.42917, -0.37727,0.80643) 

3 1 1 

74204248079595 1034415034603 307388467818' 
0 60087131 218999921 

V 0 14593 -16782 , 

j 4204248079595 (0.03212,0.32161,-0.35373) 



^3063324517380 662309451492 -144083620104 
0 132428244 106307652 

^ 0 10287 -14874 

j 3063324517380 (-0.99628,-0.39485,1.3911) 


/ 938139713086 -448140560830 27116000512 \ 

0 3374186 -90561028 

\ 0 10388 -773 / 

938139713086 (1.09384,-0.61605,-0.47780) 


o 

O 

I 
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